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We have searched for neutrino-less r lepton decays into £ and V°, where £ stands for an electron 
or muon, and V° for a vector meson (0, uj, K*° or K*°), using 543 fb _1 of data collected with the 
Belle detector at the KEKB asymmetric-energy e + e~ collider. No excess of signal events over the 
expected background is observed, and we set upper limits on the branching fractions in the range 
(0.7-1.8) x 1(T 7 at the 90% confidence level. These upper limits include the first results for £lu as 
well as new limits that are 3.6 — 9.6 times more restrictive than our previous results for £(j), £K*° 
and £K*°. 

PACS numbers: 



INTRODUCTION 

In the Standard Model (SM), lepton-flavor- violating 
(LFV) decays of charged leptons are forbidden; even 
if neutrino mixing is taken into account, they are still 
highly suppressed. However, LFV is expected to appear 



in many extensions of the SM. Some such models pre- 
dict branching fractions for r LFV decays at the level 
of 1(T 8 - 1CT 7 [E i, d, which can be reached at the 
present B-factories. Observation of LFV will then pro- 
vide evidence for new physics beyond the SM. 

In this paper, we report on a search for LFV in t~ de- 
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cays into neutrinoless final states with one charged lep- 
ton £~ and one vector meson V°: e~6, e~ui, e~K*°, 
e~K*°, [i~<f>, n~ui, [i~K*° and [i~K*° A search for 
the l~(j>, £~K*° and £~K* Q modes was performed for 
the first time at the CLEO detector, where 90% con- 
fidence level (CL) upper limits (UL) for the branching 
fractions in the range (5.1 — 7.5) x 10~ 6 were obtained 
using a data sample of 4.79 fb" 1 [5]. Later we carried 
out a search for these modes in the Belle experiment us- 
ing 158 fb _1 of data and set upper limits in the range 
(3.0 - 7.7) x 10~ 7 0. Here we present results of a new 
search based on a data sample of 543 fb _1 corresponding 
to 4.99 x 10 8 r-pairs collected with the Belle detector [7j 
at the KEKB asymmetric-energy e + e~ collider 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector, a 50- 
layer central drift chamber, an array of aerogel threshold 
Cherenkov counters, a barrel-like arrangement of time- 
of-flight scintillation counters, and an electromagnetic 
calorimeter comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside the coil is 
instrumented to detect if£ mesons and identify muons. 
The detector is described in detail elsewhere [7[. Two 
inner detector configurations were used. A 2.0 cm ra- 
dius beam-pipe and a 3-layer silicon vertex detector were 
used for the first sample of 158 fb -1 , while a 1.5 cm ra- 
dius beam-pipe, a 4-layer silicon detector and a small-cell 
inner drift chamber were used to record the remaining 
385 ib" 1 0]. 

EVENT SELECTION 

We search for t — > £<p, £lu, £K*° and £K*° candidates 
in which one r decays into a final state with a £, two 
charged hadrons (3-prong decay) , and the other r decays 
into one charged particle (1-prong decay), any number of 
7's and missing particle(s). We reconstruct cf> candidates 
from K+K~, tu from tt+tt-tt , K*° from K+n' and K*° 
from K~n + . 

The selection criteria described below are optimized 
from studies of Monte Carlo (MC) simulated events and 
the experimental data in the sideband regions of the AE 
and Mi nv distributions described later. The background 
(BG) MC samples consist of t+t~ (1524 fb -1 ) generated 
by KKMC [lOj . qq continuum, and two-photon processes. 
The signal MC events are generated assuming a phase 
space distribution for r decay. 

The transverse momentum for a charged track is re- 
quired to be larger than 0.06 GeV/c in the barrel re- 
gion (—0.6235 < cos 9 < 0.8332, where 9 is the polar 
angle relative to the direction opposite to that of the in- 
cident e + beam in the laboratory frame) and 0.1 GeV/c 
in the endcap region (—0.8660 < cos 9 < —0.6235 and 
0.8332 < cos 6* < 0.9563). The energies of photon can- 



didates are required to be larger than 0.1 GeV in both 
regions. 

To select the signal topology, we require four charged 
tracks in an event with zero net charge, and a total en- 
ergy of charged tracks and photons in the center-of-mass 
(CM) frame less than 11 GeV. We also require that the 
missing momentum in the laboratory frame be greater 
than 0.6 GeV/c, and that its direction be within the de- 
tector acceptance (—0.8660 < cos 6* < 0.9563), where the 
missing momentum is defined as the difference between 
the momentum of the initial e + e~ system, and the sum 
of the observed momentum vectors. The event is subdi- 
vided into 3-prong and 1-prong hemispheres with respect 
to the thrust axis in the CM frame. These are referred 
to as the signal and tag side, respectively. We allow at 
most two photons on the tag side to account for initial 
state radiation, while requiring at most one photon for 
the £(j), £K*°, £K*° modes, and two photons except for 7r° 
daughters for the £ui modes on the signal side to reduce 
the qq BG. 

We require that the muon likelihood ratio P M be 
greater than 0.95 for momentum greater than 1.0 GeV/c 
and the electron likelihood ratio P e be greater than 0.9 
for momentum greater than 0.5 GeV/c for the charged 
lepton-candidate track on the signal side. Here P x is the 
likelihood ratio for a charged particle of type x (x = ji, 
e, K or 7r), defined as P x = L x / (J2 X L x ), where L x is the 
likelihood for particle type x, determined from the re- 
sponses of the relevant detectors [ll| . The efficiencies for 
muon and electron identification are 92% for momenta 
larger than 1.0 GeV/c and 94% for momenta larger than 
0.5 GeV/c. 

Candidate <j) mesons are selected by requiring the in- 
variant mass of K + K~~ daughters to be in the range 
1.01 GeV/c 2 < M K+K - < 1.03 GeV/c 2 . We 
require that both kaon daughters have kaon likeli- 
hood ratios Pjc > 0.8 and electron likelihood ra- 
tios P e < 0.1 to reduce the background from e + e~ 
conversions. Candidate lu mesons are reconstructed 
from 7r + 7r~7r with the invariant mass requirement 
0.757 GeV/c 2 < M n+n -„o < 0.808 GeV/c 2 . The 
7T° candidate is selected from 7 pairs with invariant mass 
in the range, 0.11 GeV/c 2 < AL n < 0.15 GeV/c 2 . 
In order to improve the u mass resolution, the tt° 
mass is constrained to be 135 MeV/c 2 for the u 
mass reconstruction. Candidate K*° and K*° mesons 
are selected with K^n^ invariant mass in the range 
0.827 GeV/c 2 < M Krr < 0.986 GeV/c 2 , and requir- 
ing that the kaon daughter have Pk > 0.8 and both 
daughters have P e < 0.1. Figs. [T{a,b,c) show the invari- 
ant mass distributions of the (f>, u> and K*° candidates 
for t~ — » fi~(j), t~ — * fi~uj and t~ — » (i,~K , respec- 
tively. The estimated BG distributions agree with the 
data. The main BG contribution is due to qq events 
with (j) mesons for the t~ — > £~<f) mode, t~ — > n~ui> T 
with the pion misidentified as a lepton for the r~ — > £~lj 
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TABLE I: Selection criteria using p m iss(GeV/c) and 
m miss((GeV/c 2 ) 2 ) where p m iss is missing momentum and 
w-miss is missing mass squared. 



Mode 


Selection criteria 


T 


- -*■ e-<t> 


Pmiss > |m mi8S and m miss > -0.5 


T 


- -» e~u 


Pmiss > 3»Tlmiss ~ 3 an d m miss > ~ 0.5 


T~ 




Pmiss > ~4~5 m miss ~ g and p m i ss > 8m miss 


T~ 




Pmiss > 575^miss YT ^rniss ^ 


T~ 


-> fi-K*° 


Pmiss > gT^Tlmiss and 77l miss > —0.5 


T~ 




Pmiss > ^^miss and p m iss > YT4^ m i ss 



mode, and r — > 7r 7r + 7r with one pion misidcntificd 
as a kaon and another misidcntified as a lepton for the 
t~ -» £-K*° and f K*° modes. 




tracks and 7's on the signal side and that on the tag side, 
as cos a > —0.999, cos a > —0.996 and cos a > —0.990, 
respectively. 

To identify signal r decays, we reconstruct the invari- 
ant mass of £V°, M; nv , and the energy difference, AE, be- 
tween the sum of energies on the signal side and the beam 
energy -Ebcam in the CM frame. Signal events should be 
distributed around M- mv = M T and AE = 0, where M T is 
the nominal r mass. For the £lu modes, we use the beam 
energy constrained mass M\, c as the invariant mass Mi nv , 
where M^ c = y/ E^ cam — (p T ) 2 , in order to improve the 
mass resolution, which is smeared due to the 7 energy 
resolution. For the calculation of the r momentum p T , 
we replace the magnitude of the it momentum with the 
momentum calculated from the beam energy and the en- 
ergies of charged tracks on the signal side, while we fix 
the direction of the ir° momentum. 

The resolutions in AE and Mj nv evaluated using the 
signal MC are summarized in Table [TTJ We define the 
signal region in the AE — Mj nv plane as a ±3c ellipse. In 
order to avoid biases in the event selection, we blinded 
the signal region until the analysis is finalized and used 
the data in a ±10ct sideband box to estimate BG. 



TABLE II: Resolutions in Minv in MeV/c 2 and AE in MeV. 
The superscripts low and high indicate the lower and higher 
sides of the peak, respectively. 



,n nil I 

1 1.2 . 

(Gev/c z ) 





Mode 


h lg h 


low 


high 


low 


r 


- -> fj,-(j) 


3.4 ±0.2 


3.4 ±0.2 


13.2 ± 0.4 


14.0 ±0.5 


r 


— > e~ 4> 


3.7 ±0.1 


3.6 ±0.1 


13.3 ±0.7 


15.4 ±0.7 


T~ 


— » fl^UJ 


5.9 ±0.1 


6.2 ±0.1 


19.3 ±0.6 


30.3 ±0.8 


T 


— > e~to 


6.1 ±0.1 


6.5 ±0.1 


20.4 ±0.7 


32.5 ± 1.3 


T~ 


-> iTK*° 


4.5 ±0.4 


4.5 ±0.4 


13.8 ±0.3 


14.4 ±0.4 


T~ 


-» e-K*° 


4.3 ±0.1 


5.1 ±0.1 


12.9 ±0.3 


18.0 ±0.4 


T~ 




4.7 ±0.1 


4.4 ±0.1 


14.0 ±0.3 


15.0 ±0.3 


T~ 


-» e-R*° 


4.6 ±0.1 


4.9 ±0.1 


12.6 ±0.6 


17.8 ±0.5 



FIG. 1: The mass distribution of (a) 4> — > K + K~ for t~ — > 
fj,~ (f>, (b) uj — > 7r + 7r _ 7r° for r~ — > pT lo and (c) K*° — > K + ir~ 
for r~ — > fi~ K*° after muon identification. The points with 
error bars are data. The open histogram shows the expected 
t + t~ BG MC and the hatched one qq MC and two-photon 
MC. The regions between the vertical red lines are selected. 



To reduce the remaining BG from t + t~ and qq, we re- 
quire the relations between the missing momentum p m ; ss 
(GeV/c) and missing mass squared m^ iss ((GeV/c 2 ) 2 ) 
summarized in Table [U 

For the £uj (£K*° and £K*°) mode, we require that 
the opening angle between the lepton and u> (K*°) on 
the signal side in the CM frame, 6>£f (6>g^„ ), satisfy 
cosflgf < 0.88 (cos 6»™ o < .93), respectively. _ To 
remove two-photon BG for the e<j>, eoj and eK*°(K*°) 
modes, we add a condition on an opening angle, a, be- 
tween the direction of the total momentum of charged 



RESULTS 

After all selection requirements, a few events remain 
in the signal region, as shown in Figs. [2J [3J 5] and [3] 
For the £<j), £K*° and [iK*° modes, the expected number 
of BG events in the signal region are estimated using the 
sideband data, assuming that the BG distribution is flat 
in the ±10cr box. For the fiu; modes, we estimate the 
BG contribution in the signal region using the BG MC 
distribution normalized to the ratio of data and MC in 
the sideband region, because a rather large number of BG 
events remain, which are mainly from r~ — * tt~lov t . For 
the t~ — > e~tjj and e~K*° modes, since no events remain 
in the ±10(7 box for data, the number of expected BG 
events in the signal region is zero. We calculate the error 
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FIG. 2: AE vs. Minv distributions for (a) r~ — > jj,~(f) and (b) 
r~ — > e~ <j>, after all selection criteria. Dots are data, yellow 
boxes show the signal MC, purple dots are background from 
r — > 4>np MC and open circles show other r-pair backgrounds. 
The elliptical area is the 3<r signal region. 
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FIG. 5: AE vs. M nv distributions for (a) r -> /x A"* and 
(b) r~ — > e~A*°, after all selection criteria. Dots are data, 
yellow boxes show the signal MC and open circles show r-pair 
background MC. The elliptical area is the 3<r signal region. 
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FIG. 3: AE vs. Mmv distributions for (a) r~ — » [m~lv and 
(b) r~ — » e~oj, after all selection criteria. Dots are data, 
yellow boxes show the signal MC, and open circles show r-pair 
background MC. The elliptical area is the 3<r signal region. 



TABLE III: Number of events in the sideband region. 



Mode 


A^v in the ±10a box 
in the sideband region 
data MC 


T~ 


- -» ix 4> 


2 


1.68 ± 1.17 


T~ 


— ► e~(j) 


2 





T~ 


— > /J,~U! 


8 


8.97 ± 1.72 


T~ 


—* e~u> 





1.07 ±0.62 


T~ 


-» 


4 


3.94 ± 1.61 


T~ 




1 


1.58 ± 1.11 


T~ 


-» 


2 


0.90 ± 1.11 


T~ 







0.06 ± 0.06 



from the number of remaining MC events in the ±10<7 
box assuming a flat BG distribution. The number of 
events in the ±10ct box excluding the signal region and 
the number of the expected BG events are listed in Tables 
IIIII and HVl The comparison between the data and MC 
shows reasonable agreement; the BG is suppressed well 
by the event selection. 
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0.05 


-0.05 
-0.1 



1 0.1 
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-0.05 
-0.1 
-0.15 



Ibl 
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M uK .(GeV/c 2 ) 
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M eK .(GeV/c 2 ) 



FIG. 4: AE vs. Minv distributions for (a) r~ -> (i~ K*° and 
(b) r~ — > e~ K*° , after all selection criteria. Dots are data, 
yellow boxes show the signal MC and open circles show r-pair 
background MC. The elliptical area is the 3<r signal region. 



From the remaining number of data events in the sig- 
nal region and the number of expected BG events, we 
evaluated the upper limit sgo on the number of signal 
events at 90% CL with systematic uncertainties included 
in the Feldman-Cousins method [13] using the POLE 
code 13(. The main systematic uncertainties on the de- 



tection efficiency come from track reconstruction (1.0% 
per track), electron identification (2.2%), muon identifi- 
cation (2.0%), kaon identification (1.4% for (f> reconstruc- 



TABLE IV: Number of expected BG events. 



Mode 


A ev of expected BG 
in the signal region 


r 




0.11 ±0.08 


T 


— > e~ (j> 


0.11 ±0.08 


r" 




0.19 ±0.21 


r 


— ► e~to 


0±0.07 


T~ 


-> n~K*° 


0.22 ±0.11 


r~ 


-> e-R*° 


0.05 ±0.05 


r~ 


-» fi-K*° 


0.10 ±0.07 


T~ 


-» e-R*° 


0±0.02 
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tion and 1.1% for K* ), ir° reconstruction (4.0%), statis- 
tics of the signal_MC (1.3% for £(f>, 0.7% for £co and 0.6% 
for £K *° and £K*°) and uncertainties in the branching 
fractions for 4> — > K + K~ and uj — > 7r + 7r~7r° (1.2% and 
0.8%). The uncertainty in the number of r-pair events 
mainly comes from the luminosity measurement (1.4%). 
Upper limits on the branching fractions B are calcu- 



lated as B < 



2N TT <i- 



where N TT = 4.99 x 10 s , which 



we calculate using cross section of 0.919 nb according to 
141 ] . is the total number of the r-pairs produced and e 
is the signal efficiency including the branching fractions 
of -> K+K~, uj -> tt+tt-tt and -> K+ir- [H. 
The resulting upper limits on the branching fractions are 
summarized in Table fVl 



TABLE V: Summary of the number of observed events iV bs, 
detection efficiency e, total systematic error Ae/e, 90% CL 
upper limit of signal events sgo and 90% CL upper limit of 
branching fractions. 





Mode 


iV ob8 


e 


Ae/e 


S90 


UL on BF 








(%) 


(%) 




(90% CL) 


t" 


" -»■ M~<^> 


1 


3.14 


5.2 


4.21 


1.3 x 10" 7 


r" 


—> e~ ' (j> 





3.10 


5.3 


2.34 


7.6 x 10~ 8 


T~ 







2.51 


6.3 


2.25 


9.0 x 10~ 8 


T~ 


— » 


1 


2.46 


6.3 


4.34 


1.8 x 10~ 7 


T~ 


-> /x^*° 





3.71 


4.8 


2.24 


6.1 x 10" 8 


T~ 







3.04 


4.9 


2.42 


8.0 x 10~ 8 


T~ 


-> /x^*° 


1 


4.02 


4.8 


4.23 


1.1 x 10~ 7 


T~ 







3.21 


4.9 


2.45 


7.7 x 10" 8 



SUMMARY 

We have searched for LFV decays r~ — > £~<f), £~uj, 
£~K*° and £~K*° using a 543 fb -1 data sample from 
the Belle experiment. No evidence for a signal is observed 
and upper limits on the branching fractions are set in the 
range (0.6-1.8) x 10~ 7 at the 90% confidence level. This 
analysis is the first search for t — > £u> modes. The results 
for the t~ -> £-(/), £~K*° and f K*° modes are 3.6-9.6 
times more restrictive than our previous results obtained 
using 158 fb -1 of data. The sensitivity improvement in- 
cludes a factor of 3.4 in data statistics and an optimized 
analysis with higher efficiency and much improved BG 
suppression. The improved upper limits can be used to 
constrain the parameter spaces of various scenarios be- 
yond the SM. 
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